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We study cosmological perturbations in the brane models with an induced Einstein-Hilbert term 
on a brane. We consider an infiaton confined to a de Sitter brane in a five-dimensional Minkowski 
spacetime. Inflaton fluctuations excite Kaluza-Klein modes of bulk metric perturbations with mass 
= -2(21 - 1)(^ + 1)H'^ and = -2^(2£ + 3)//^ where I is an integer. There are two branches 
(± branches) of solutions for the background spacetime. In the + branch, which includes the self- 
accelerating universe, a resonance appears for a mode with nn? = 2H^ due to a spin-0 perturbation 
with = 2H^ . The self-accelerating universe has a distinct feature because there is also a helicity- 
mode of spin-2 perturbations with = 2H^ . In the — branch, which can be thought as the 
\^ • Randall-Sundrum type brane-world with the high energy quantum corrections, there is no resonance. 

' At high energies, we analytically conflrm that four-dimensional Einstein gravity is recovered, which 

, is related to the disappearance of van Dam-Veltman-Zakharov discontinuity in de Sitter spacetime. 

■ On sufficiently small scales, we confirm that the lineariaed gravity on the brane is well described 

by the Brans-Dicke theory with ui = SHvc in — branch and ui = —SHvc in + branch, respectively, 
which confirms the existence of the ghost in -I- branch. We also study large scale perturbations. In 
+ branch, the resonance induces a non-trivial anisotropic stress on the brane via the projection of 
Weyl tensor in the bulk, but no instability is shown to exist on the brane. 
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I. INTRODUCTION 



> 

' There has been tremendous interest over the last several years in the brane-world scenario where_we are assumed 
] to be living on a four-dimensional hypersurface (brane) in a higher-dimensional spacetime (bulk) l]. The simplest 
example of this scenario is proposed by Randall and Sundrum where there is a brane embedded in a five-dimensional 
anti-de Sitter (AdS) spacetime They showed that the four-dimensional gravity is recovered at low energies, because 
gravity is confined to a single positive-tension brane even if the extra dimension is not compact. Therefore, in order 
to see the deviation from the conventional four-dimensional gravity, it is necessary to investigate the phenomena 
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I , occurred at high energies. 

' One of the most promising tools to extract the information of the extra-dimension is the primordial density fluc- 
tuations generated in the period of inflation in the early universe. In the inflation model where inflation is driven 
by an inflaton field confined to the brane, the amplitude of the curvature perturbation is calculated in the extremely 
slow- roll limit where the coupling between inflaton field fluctuation and bulk perturbations can be neglected Q . This 
^ work is proceeded to go beyond the zero-th order slow-roll approximation by solving the bulk metric perturbations 
' classically ^ |^. It is shown that we could have significant effects from the back-reaction due to the coupling to 
five-dimensional gravitational perturbations. On the other hand, it was pointed out that the localized matter on a 
brane can induce gravity on the brane via quantum loop corrections at high energies This induced gravity can 
act as the ultra-violet cut-off for the inflaton perturbations. 

Based on the induced gravity scenario, Dvali, Gabadadze, and Porrati (DGP) proposed a brane-world model with 
induced gravity in which the four-dimensional brane is embedded in a five-dimensional Minkowski spacetime. (For 
a review of the phenomenology of DGP model, see 1^.) In the presence of the induced gravity, according to the 
embedding of the brane in the bulk, there appear two branches of background solutions. In -I- branch, the size of the 
four-dimensional hypersurface becomes minimum at the brane while, in — branch, it becomes maximum at the brane. 
Inflation models in this scenario were studied in |^, 0, ^| . 

In this paper, we study the behavior of five-dimensional metric perturbations excited by the inflaton perturbations 
confined to the brane in DGP model. The solutions in — branch are regarded as the high energy limit of the Randall- 
Sundrum model including the curvature corrections on a brane, depending on the parameters. Then our analysis 
can be used to discuss the generation of primordial fiuctuations in these models. The solutions in -I- branch are 
largely used to explain the present cosmic acceleration without introducing dark energy in the DGP model (l^], but 
it turns out that they are plagued by the ghost instability 13, IJ^ J5j_16, 17]. However, the study of the gravitational 
property of this model in the presence of matter and cosmic expansion is still very limited 0, ITqI l20j | so it is also 
worth considering the solutions in this branch. 

The structure of the rest of the paper is as follows. In section II, we derive the basic equations for the five- 
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dimensional metric perturbation and the boundary condition imposed at the brane and then summarize the solution 
for the vacuum brane obtained in [T5l |. In section III, we provide the solutions for curvature perturbations on the 
brane in the presence of the scalar field on the brane. Next we discuss the appearance and disappearance of van Dam- 
Veltman-Zakharov discontinuity in section IV. In section V we show the comparison of the behavior of perturbations 
on small scales with the Brans-Dicke theory. The solutions for perturbations on large scales are discussed in section 
VI. Section VII is devoted to conclusions. Technical details of the calculation of the m? = 2H^ mode and the brief 
summary of the cosmological perturbation in four-dimensional Brans-Dicke gravity are presented in the Appendix A 
and B, respectively. 



II. BULK GRAVITONS WITH A DE SITTER BRANE 



A. Background cosmology 

We consider a four-dimensional brane- world model with a five-dimensional infinite- volume bulk. The action is given 
by 



if + I/branc + „ 5" 



(1) 



where is the fundamental five-dimensional gravitational constant, K is the trace of extrinsic curvature Kab on 
the brane and Lbianc is a Lagrangian for brane-localized matter. In addition, we consider the last term, an intrinsic 
curvature term on the brane which plays a crucial role in this model. Here, k| = SttG is the four-dimensional 
gravitational constant. We define a crossover scale Tc as 

We are interested in the gravitational property of this model in a cosmological setting. By assuming the flat 
Friedmann-Robertson- Walker metric on the brane and neglecting the effect from the five-dimensional Weyl tensor in 
the background, the following Friedmann like equation is derived 

H^--H= f p, (3) 

where iJ is a Hubble parameter on the brane and e = ±1, which is the parameter related to the embedding of the 
brane in the bulk. We call the case with e = 1 -I- branch, while t — — branch. 

In this paper, we restrict ourselves to de Sitter background H = constant. In this case, the energy density is 
constant and given by 

(4) 

The 5D solution for the metric with the four-dimensional de Sitter brane can be obtained as 

ds" - dy^ + iV(y)27^,c^x'^da;^ N[y) = 1 + ei/y, (5) 

where 7^1^ is the metric for the de Sitter spacetime and the brane is located at ?; = 0. During a slow-roll inflation, the 
Hubble scale is varying slowly, thus we can approximate the background geometry by |(SJ|. 



B. Master variable for perturbations in the Minkowski bulk 



Now let us consider the scalar perturbation. In order to solve the perturbations in this background, it is convenient 
to use five-dimensional longitudinal gauge, given by 



ds^ = (1 + 2Ayy)dy^ + 2N{y)Aydydt + N{yf - (1 + 2A)dt^ + a{tf{l + 2TZ)Sijdx'dx^ 



(6) 



where a(t) = exp{Ht). 
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In the absence of bulk matter perturbations, five-dimensional perturbed Einstein equations ^^^SG^^ = are solved 
in a Minkowski background if the metric perturbations are derived from a 'master variable', Vl (see also [2^'): 
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(7) 



where dot and prime denote derivatives with respect to t and y, respectively. The perturbed five-dimensional Einstein 
equations yield a single wave equation governing the evolution of the master variable Q in the bulk: 



- , N' 

n - 3Hn - in" - 2 — n' 

N 



,n = 0. 



(8) 



Solutions for the master equation can be separated into eigenmodes of the time-dependent equation on the brane and 
bulk mode equation: 



(9) 



where 



3Hg„ 



9m = 0, 



N' 77? 2 

f" - 2— f + — f - 



(10) 

(11) 



C. Boundary condition for Q, 



In order to solve fi, we must specify the boundary condition for Vl from the junction conditions at the brane. The 
junction conditions in this model have been found already in literatures (see |23,|23), but, we re-derive these so that 
the property of linearized gravity is transparent. 

For simplicity, we consider the case where the anisotropic stress of the matter on the brane does not exist. There 
are two important contributions to the junction conditions. One is that, from Eq. H83|) . we must take into account not 
only the matter energy momentum, but also the contribution from the induced gravity on the brane. The other is that 
we should take into account the brane bending ^ |25l l26j | . The metric perturbations in five-dimensional longitudinal 
gauge evaluated at y = is not the induced metric perturbations on the brane. This is due to the anisotropic stress 
coming from the induced gravity. The junction conditions are given by 



T + 2HT 



1 
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- 2e - AHi + 6H^^ - TT^^ = i^^Spn, + Sp, 



3a2 



--T + 2i- 2HC = ^i'iSqm + Sq^), 



7,2 ■ 



mi - m'^i + 2—i = n'idp^ + 5p^), 



(12) 
(13) 
(14) 



where 



T = n' - eHn. 



(15) 



Here the brane bending mode is given by 



C = -rc(A(b) -l-7^(fc)), 



(16) 
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where and 7?.(t,) are the induced metric perturbations on the brane that are related with the metric perturbations 
in five-dimensional longitudinal gauge as 



1 



n 



1 - 2eHrc 
1 



(6) 



1 - 2eHrr 



{l-eHrc)A + eHrcn 
(1 - eHra)n + eHr^A 



(17) 
(18) 



where TZ and A are metric perturbations in five-dimensional longitudinal gauge evaluated at y = 0, which can be 
expressed by O as 



ba 



3eHT - 3H{n - HQ) + — O 



A = 



1 r 

6a 



- 3eHT - 30 + GHil - 3H^n - 2— ft 



(19) 
(20) 



dpm-, Spm and Sqm are the components of the perturbed energy momentum tensor of the matter on the brane and 



5q„ ^ 4 



-3i/2A(fc) - + ^(6) + 3i^7^(b) + ^(^(fc) + n^b)) 



Spg = 



(-ff^(b) - 7?.(fc))i 
3i^(-7^(b) + i^A(fc))-— 7e(b) 



(21) 
(22) 
(23) 



are contributions from the induced gravity. If we do not take into account the induce gravity term which corresponds to 
the limit Vc oo, the quantities ^, Spg, Sqg and Spg can be neglected and the junction conditions in Randall-Sundrum 
type brane- world 0,123, 123 are reproduced. 
Defining Q as 



g = {1- 2eHrc)T - rc{2H^ - m^)n, 
the junction conditions ifT^ . ((T^ . ((TH) are amazingly simplified as 

fc2 



-3Hg 



-g = K aSp,^, 



-g = n'^aSq^, 
g + 2Hg = K^aSpn 



(24) 

(25) 

(26) 
(27) 



D. Solutions for a vacuum brane 



Before discussing the case with the matter perturbations on the brane, we briefiy summarize the case for a vacuum 
brane, ST^^, = considered in |l5l |l6l| . We will only consider a positive tension brane which means Hrc > 1 for -|- 
branch and Hrc > for — branch. The solutions obtained here will serve as homogeneous solutions when we include 
matter perturbations. 

For the vacuum brane, the boundary conditions for ft given by Eqs. (|25|l . I|26|) and H27|) reduce to a single boundary 
condition on the master variable, 

g = 0. (28) 

For both branches, we find a tower of continuous massive modes starting from — (9/4)_ff^, which are the Kaluza- 
Klein (KK) modes of the spin-2 perturbations. 

In -|- branch (e — 1), in addition to continuous modes, there are two discrete modes which are normalizablc. One is 
the mode with — 2H^ which is the spin-0 perturbation, and the other is the mode with — H^{3Hrc—l){Hrc)^^ 
which is the helicity-0 excitation of the spin-2 perturbations |2^. For Hrc > 1, the mass is in the range 
< TO^ < 2H^. If Hrc — 1, the helicity-0 excitation of the spin-2 perturbation has mass = 2H^, which 
coincides with the mass of the spin-0 perturbation. Then, there is a resonance between these two modes. It was also 
shown that the spin-2 helicity-0 mode becomes a ghost for Hrc > 1. In Figure 1, we summarize the mass spectrum 
for 1 < Hrc in + branch. 
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FIG. 1: Summary of the mass spectrum of the scalar perturbations in + branch H^ . Spin-2 perturbation has continuous modes 
with m? > (9/4)_ff^ and a discrete mode with = while spin-0 perturbation has = 2H^. In the hmit Htc — > 1, both 
the hehcity-0 excitation of spin-2 perturbation and the spin-0 perturbation have mass = 2H^ and there is a resonance [l^. 



III. SCALAR FIELD ON THE BRANE 



A. Bulk scalar modes 



Until now, we have not specified the matter on the brane. In the following, we model the matter on the brane as 
a canonical scalar field with potential V{(t>) whose homogeneous part gives approximately de Sitter universe on the 
brane jUl]- 

In this case, the junction conditions H12I) . H13|) . H14|l reduce to 

1.2 . . 

-3Hg -^g^ K^a[(l)S(j) + V,M], (29) 

g = K^a^Scj), (30) 
g + 2Hg ^ K^a[^6<l)~V,^S(j)]. (31) 

These equations can be thought as the boundary conditions for il. Combining the junction conditions (|29|l . H30|l . 
(|ST|l . we get an evolution equation for g-. 

g- (^H + 2^g + ^g = 0. (32) 

This gives the boundary condition for the time dependence of the master variable il. Assuming that (j) is slow-rolling, 
so that \(t)/4>\ H in Eq. H32|) . which is valid for the de Sitter universe, the solution for g is 

~kr] —krj 
We use the formulae for summation of Bcssel functions, 

E(-l)^f2.+ i).-ij,,,.(.) . (34) 



These show that an infinite sum of mode functions 



9 Tfi^ 

9m = {~'kr]y^''^J^{-kr]), where i/^ ^ _ _ (35) 
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can satisfy the boundary condition imposed on 5, where the spectrum of KK modes is given by 

9 

Tfl 

— = -2(2£-l)(^+l) for Ci, 



m 



-2l{2l + i) for Ca. 



(36) 



Figure 2 shows the mass spectrum of the scalar perturbation supported by the scalar field on the brane. Since there 
is a mode with m? — 2H^ regardless of the value of Hrc in the + branch, the resonance inevitably appears. This 
gives qualitative difference between — branch and + branch. In addition, for Hrc — 1, that is, on the self-accelerating 
universe, the spin-0 and spin-2 modes are already degenerates. Thus the behavior of gravity is again qualitatively 
different from Hrr > 1. 



Spectrum 



Ci 




4 




< 


► 




• 
• 
• 















spin 










4 


< 








-lOH^ 



FIG. 2: Summary of the mass spectrum of the scalar perturbation sourced by the scalar field on the brane. For the modes 
Ci, there is a mode with — 2H^ regardless of the value of Hrc- Since there is another mode with — 2H^ by the spin-0 
perturbation in the + branch, there is a resonance. In addition, for Hvc ~ 1, the helicity-0 mode of spin-2 perturbations also 
has = 2H^. 



In order to examine the behavior of the master variable Q, we will consider — branch (e = —1) and + branch (e = 1) 
separately in the following. 

— branch (e = 1) 

Since there is a horizon at y — 1/H in this branch, we should choose the solution in the y-direction so that the 
metric perturbations remain small as y 1/H. By imposing the normalizability condition on /,„, 



l/H 

Niy)-^\fmiy)\^dy < oo, 
we obtain the solution for O in the bulk (for y > 0) as 



with N{y) = 1- Hy, 



H 



21- 



(1 - Hyf+^' 



H 



1=0 



2 J {2£+l){2{£+l)Hrc + l} 
(1 - Hyf+^' 



2 J 2{i+l){{2i + 3)Hrc + l} 



(37) 



(38) 
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+ branch (e = 1) 

In + branch, the bulk is infinite. We should choose the solution in the y-direction so that the metric perturbations 
remain small as y — > cx3 and the mode is normalizablc for a single brane. We impose the normalizability condition on 
fm as 

/•CO 

N{yr^\f,n{y)fdy<^, with N{y) = l + Hy. (39) 

Because of the existence of the resonance, we must treat the mode with = 2H^ separately. We also treat Hvc = 1 
case separately as the spin-2 and spin-0 degenerate (see Appendix A for derivations). 
We obtain the solution for in the bulk (for y > 0) as 

H f-^^ ' \ 2j 2£{{2e-l)Hrc + l}^ ' ' ' 



where 



H ' \ 2 J {2£+l){2mrc + l) 

(40) 



^(m^=2m) = '^^ 2H(Hr^^- l) ^^^^^ ^{a{-kT])fim{-kij) + /3{-kij)cos{-kij)} 

^ {1 + Hy) + Hy) . , ,_2 • / , n 
-Ci hIht^-I) ^ -^"^^ 'sm(-M, (41) 

d{-kfj){-kfjy^ sin{-2kfj), (42) 
Pi-k-q) = / d(-fc^)(-fc77)-2{cos(-2fc^) - 1}, (43) 



for Hvc 7^ 1 and 



C 

^{m^=2H^) = -^(1 + Hy){-kTj)~'^{&{~kr]) sm{~kr]) + f3{-h]) cos(-fc?7)} 
C 

H -{1 + Hy) ln(l + Hy){-kT]y^{a{-kri) sm(-kri) + /3{-kri) cos{~kri)} 

2H 

-^(1 + Hy){ln{l + Hy)Y{-ki^)-^ sin(-fcr;), (44) 

(45) 
(46) 
(47) 
(48) 
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— (3{—kff) sin(— 2fc?7)], 



for Hrr = 1. 



B. Curvature perturbation on the brane (general solutions) 

Since we have obtained the general solutions for 17, from Eqs. (I17|) . H18|l . (|19|l and (|2()|l the induced metric pertur- 
bations on the brane can be obtained. We concentrate on the behavior of the curvature perturbation on the brane 
7?,(f,). As in the previous section, we will consider — branch (e = —1) and + branch (e = 1) separately. 
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— branch (e = —1) 

In — branch (e = —1), the curvature perturbation on the brane is given as 
\/2^Ci ^ (-l)^(2£+i) 



'^^^^ 6 k i2Hr,, + l)^^i2£+l){2i£+l)Hr^ + l} 



x[6(^+l)(2^+l)i/re(-A;r7)"'-/2£+i(-M + 2(2£+l)(-fcr7)'^2^-i(-M-(-fc^)^'/2£-i(-fc^)] 
^ (-1)^(2^+1) 



6 k {2Hr^ + 1) ^ + l){{2t + i)Hr^ + 1} 
X W + l){2l + 3)Hr,{-krj)--^J,,^s{~kfj) + 2{2£ + 2){~krj)-^ J,,^i{-kTj) - (-ki^)i J,,^i{-kTj)]. 

(49) 

+ branch (e = 1) 

In + branch (e — 1), the curvature perturbation on the brane is given as 

V2^Ci (-1)^(2^+ ^) 

6 k {2Hr, - 1) ^ 2i{(2l - l)Hr, + 1} 
X [{U{21 - l)Hr, + 12^ + 3}(-fcr;)-i J2,+i(-fcr;) - 2(2^ + l){-ki^)-- J^,_^_{~kl^) + {-kr,)i J^,^-^{-ki^)] 
_^ (-l)'(2^+§) 



6 /fc (2i?rc- 1) ^ (2^+l){2^i/rc + l} 



where 



x[{6^(2^+ l)iJr, + 12^ + 9}(-fcr;)-^ J2,+|(-fcry) - 4(€+ l)(-M'^2^+i(-fc'7) + (-fc'7)^^2£-i(-fc'?)], 

(50) 



7e(.)(™.=2^.) = Ty(^^[3(ffr.-l)(-fc,)-isin(-fc,) 



2 



{—3{—kT]) sin(— fcr;) + 3 cos{—kri) + {—kr]) sin(— fc??)} 
{— 3(— fc77)^"^ cos(— fcr;) — 3 sin(— fcry) + {—kri) cos(— fcry)}] , (51) 



for Hvc 1 and 



^(6)(m2=2ff2) = '^I2fc [~"^(^^^) ^sin(-fc?7) 

+a{—kri){—3{—kri)^^ sin(— fc?7) + 3 cos(— fcr;) + {—krj) sin(— fcry)} 
+f3{-kri){~3{-'kriy^ cos(-fcr;) - 3 sin(-fc?7) + {-kr]) cos(-fc?7)}] , (52) 

for Hr^ = l. 

We should note that there are homogeneous solutions that satisfy ^ = 0. These homogeneous solutions also induce 
curvature perturbations on the brane. 



IV. VAN DAM-VELTMAN-ZAKHAROV DISCONTINUITY 
A. Curvature perturbation on the brane in high energy limit 



Since the background cosmology given by Eq. Q recovers the conventional four-dimensional cosmology based on 
the general relativity at high energies, we expect that the gravity behaves like four-dimensional general relativity in 
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this limit. Here, we show this fact analytically using our solutions. If we take the limit Htc — > oo for fixed —krj, by 
using the formula Eq. (|34|1 . the sum of the infinite ladder of the discrete modes can be performed and we obtain 

7^(^,) = = (Ci cos(-/c7/) + C2 sin(-fc77)) , (53) 

in both branches, which agree with the solutions in general relativity Eqs. HB18|I . HB19|I . Therefore, we see that 
four-dimensional general relativity is recovered in this high energy limit. This confirms the result first obtained in 
Ref. [2|. 

In this limit, the brane bending vanishes ^ = — 7'c(A(ft) +7?.(fc)) — 0. We can understand this result in a very simple 
way as follows. In terms of A and TZ, ^(t) and 'R-(b) are given by 

-4w+7^W- l_2,g,j -4 + ^)- (54) 

Then at high energies Hr^ ^ 1, the brane bending mode is shielded by the curvature of the brane [l8l |. 

However, this argument holds only for Htc — > 00. For a finite Hr^^ if we do not fix ^krj and take large enough 
—krj, we can no longer perform the summation of the infinite ladder of the mode and obtain the solutions (|53l) . Thus 
the theory deviates from general relativity. This is related with the well-known van Dam-Veltman-Zakharov (vDVZ) 
discontinuity first discussed in the context of the Pauli-Fierz theory for a massive spin two field [^D] and also in 
DGP model 

B. van Dam-Veltman-Zakharov discontinuity 

Let us consider small scales limit —k-q ^ 1. In this region, (0,0) component of the junction conditions Eqs. 114(1 
and the solution for the brane bending give 

On small scales, we can assume Hil' <C {k'^ /a^)Vl, Then, Eq. H55I) and (|56|l becomes completely the same as the one 
obtained from Brans-Dicke (BD) theory l|B12|) and l|B15p 19], provided that 

5^ = {\lK^)t (57) 

Thus in general we do not recover four-dimensional general relativity but BD theory. This is due to the fact that 
four-dimensional gravity is recovered by a continuum of massive states of five-dimensional graviton. 

In the high-energy limit, we have essentially shown that the BD parameter becomes infinity for Hrc — + 00. This re- 
sults remind us that in the case of massive gravity, if we introduce a cosmological constant A(4), the vDVZ discontinuity 
disappears if the limits m — > is taken for fixed A(4) [sil Is^ l33l | . 

For a finite Htc, the theory is not described exactly by general relativity and this difference shows up at —/cry ^ 00. 
Then we need to know what kind of BD theory is realized for a fixed Hrc- The value of the BD parameter oj is 
obtained only by solving the equation of the motion for the BD scalar field l|Blip because it is the only equation 
which includes the BD parameter. In the BD theory, the BD scalar obeys a simple four-dimensional equation of 
motion. However, in the brane-world model, the behavior of the brane bending mode can not be obtained unless we 
solve the five-dimensional perturbations. Now, since we have a solution for fi, we can find the solutions for the brane 
bending mode. In the next section, we identify the BD parameter using the solution for the brane bending. 



V. COMPARISON WITH BRANS-DICKE THEORY 



In this section, we derive the solutions for the brane bending mode and identify the BD parameter for the effective 
theory on small scales. 
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A. Brane bending mode 

Here, we show the solutions for the brane bending mode ^ in — branch and + branch separately. 
— branch (e = — 1) 

In — branch (e = —1), the brane bending ^ is given as 



^ V2^C,H Hr, ^ (-1)^(2^+1) 
^ 3 k (2ifre + l)£;'(2^+l){2(f+l)ifre + l} 



x[3{£+l){2l+l){-kv)-^J^,^.{-kv)-2{2i+l){-kv)ij^,_i{-k^^ 

V2^C2H Hr, y. (-1)^(2£+|) 
3 k~ {2Hrc + 1) ^ + 1){{2£ + 3)Hrc + 1} 

x[3(^+l)(2^ + 3)(-fcr?)-5J2,+2(-M-2(2^ + 2)(-M^-/2^+i(-M + (-M^^2^-i(-M]- (58) 
+ branch (e = 1) 

In + branch (e = 1), the brane bending ^ is given as 

V2^CiH H ^ (-1)^(2£+1) 
3 k 2Hrc - 1 ^ 2l{{2£ - l)Hrc + 1} 

X [-3{£ + l){2e + l)(-M"^^2^+i + 2(2£ + l){-kv)ij2e-i - {-kv)^ J^t-^] 
V2^C2H Hr, ^ (-1)^(2£+|) 



3 k 2Hrc-l^^i2e + l){2£Hrc + l} 

x[-3(£+l)(2^ + 3)(-M"^^2€+i+2(2£ + 2)(-fc7?)*^2€+i-(-M^^2^-i]. (59) 

where for Hvc ^ 1 

+a(— fcry){3(— fcr/)~^ sin(— fcr/) — 3cos(— A;ry) — {—kr])sm{—kr])} 
+/3(-A;r?){3(-fc7?)"^ cos(-fcr?) + 3 sin(-A;r?) - (-fcr?) cos(-fc7?)}] , (60) 

CiH r 

^(m2=2H2) = a{- kf]) kT]y^ sm{- krj) ~ 3cos{- kij) ~ {-kr]) ami- krj)} 

+P{—kr]){3{—kri)~^ cos(— fcr?) + 3 sin(— fciy) — {—krj) coa{—kri)} 

3 3 
+-a{—k'r]){—kr])~^ sin(— /cry) + -(3{—krj){—kr])~^ cos(— /cry) 

+3(-/cr?)-^ sin(-fc7?) , (61) 

There are homogeneous solutions that satisfy Q = Q. These arc given by 

VI oc {-k'q)-^/'^Z^{-kr]), (62) 

where Z^, is a linear combination of Bessel functions and u is determined by the spectrum of massive modes. On small 
scales, regardless of the value of u, the solution for the brane bending from the homogeneous solutions are given by 

^ = di{—kr]) sm{— krj) + d2{—kr]) cos(— /cry), (63) 

where di and d2 are arbitrary constants. 



and for Hrr = 1 
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B. Perturbation of Brans-Dicke scalar 

Now we examine the behavior of in the BD theory in the presence of the scalar field 0, whose dynamics is 
completely the same as the brane-world model. Since we are interested in the behavior of b^p on sub-horizon scale, 
we take the small scale limit. Then, from Eqs. (|3UII and can be evaluated as 

(t)'d(t)' = 2 [Ci cos(-fcr;) + C2 sin(-A:ry)] , (64) 

where in order to obtain Eq. (|64|) . we used the fact that {(f>/<j)) ^ -ff <C fc/a in the de Sitter background and subhorizon 
scale. 

This term serves as a source term in the equation of motion of the perturbed Brans-Dicke scalar field IjBlip which 
can be simplified in the de Sitter background as 

Sip" + ^Sip' - V'^Sip = ~2C'^ (t>'S(j)', (65) 

where ^ is a constant related with the BD parameter a; as = (6 + Aoj)^^. Then by solving Eq. H65|) and keeping 
leading terms, we obtain 

Sip = di(—kri) siii{~kri) + d2{—k'ri) cosl—kr]) 

+ ^ [-C2{~'k?]) sin{-kf]) + Ci{-kr]) cos{~h])] Si{~2kf]) 

+ \ [-Ci{-kTj)sm{~kT]) - C2{~krj)cos{~krj)\ Ci{-2krf) 

+ - \Ci{~krj) hi{—krj) sin(— fcry) — C2{—kr]) ln(— fcyy) cos(— fcry)] , (66) 

where di and d2 are integration constants related with the homogeneous solutions and Ci and C2 are given by 

Cr = C2 = (67) 

Note that the homogeneous solutions behave in the same ways as the brane bending mode coming from homogeneous 
solutions for f2 (|63|l . These homogeneous solutions are fixed by appropriate boundary conditions. 

C. Numerical results 

Here, we identify the value of the Brans-Dicke parameter cu. For this purpose, we evaluate Eqs. H58|l . (|59|l and 
(|66|l numerically. In practice, we must approximate the infinite sum to proceed the calculations. In our previous 
paper , we checked that we can approximate the infinite summation by introducing a cut-off £c into the summation 
because Bessel function satisfies Ju{z) — > for v > z. Therefore, as long as we start from a finite time —krji, we can 
approximate the infinite ladder of the modes by introducing sufficiently large cut-off £c and the result is insensitive 
to 4. 

Here, we compare the part of the numerical solutions including Ci of Eqs. (|59|l with (|66|l . The same result applies 
to the solution including C2. For the numerical calculation, we construct the following dimensionless quantities; 

k 1 Ik 

^ = -^^(^BDl +7^BDl), (68) 

where the subscript 1 denotes that they correspond to the part including Ci. 

Figure 3 shows ^1 (solid line) and S<pi (dotted line) as functions of —krj in the scale smaller than the Hubble scale 
{—krj = 1). The Brans-Dicke parameter is chosen as 3Hrc. Here we take Hvc — 2 and choose the amplitudes of the 
homogeneous solutions for ^i^i, di and ^2 so that it agrees to ^1 with di ^ d2 — at —krj = 200 . We see that the 
behavior of perturbations on scales much smaller than the Hubble scale is well described by the Brans-Dicke theory 
with uj ~ 3Hrc- It can also be seen that, at least for — fcry < 120, the deviation of ^1 from Sipi become large and the 
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FIG. 3: ^1 (solid line) and 5(fi (dotted line) as functions of a; = —kr/. The Brans-Dicke parameter is SHr^- Here we take 
Hvc = 2 . The dotted line completely coincides with the solid line on small scales (large x). The deviation become obvious 
around x < 120. 



solution is no longer well described by the Brans-Dicke theory with u; = 3Hrc- Figure 4 also shows (solid line) and 
S(fi (dotted line) with ut = 3Hrc as functions of ~h] for Hvc — 0.1. We also choose the homogeneous solutions so 
that both coincide at —k-q = 200 in the same way. Again we see the solution is well described by the BD theory with 
Lv = ZHtc, but the deviation of S^i5 from (pi become large for ^krj < 180. This is natural because, for Hvc = 0.1, Vc 
becomes smaller than the horizon scale (corresponding to —hi] = 10). 




FIG. 4: ^1 (solid line) and 5(fi (dotted line) as functions of a: = —kri. The Brans-Dicke parameter is 3Hrc- Here we take 
Hrc = 0.1 . The dotted line completely coincides with the solid line on small scales (large x). The deviation become obvious 
around x ~ 180. 



These results are consistent with Ref. 0, . In Ref. ^| , it is shown that the linear growth of density contrast 
in non-relativistic spherical collapse of dust is described the Brans-Dicke theory with uj = SHvc for l/H 3> r when 
Hrc is greater than 1. 

Even though we do not show the results explicitly here, we can confirm that, in the + branch, similar results can be 
obtained that the perturbations can be described by the Brans-Dicke theory with w = —ZHrc on small scale. This is 
also consistent with 0, 0| . In the Brans-Dicke theory, the Brans-Dicke scalar becomes a ghost if the Brans-Dicke 



13 



parameter satisfies w + 3/2 < 0. This result implies that the brane bending ^ is a ghost for Htc > 1 in + branch. 



VI. CURVATURE PERTURBATION ON THE BRANE ON LARGE SCALES 



A. Solutions on Large scales 

In this section, we study the behavior of perturbations on large scales {—krj) <^ I. In this limit, the solution for 
is dominated by the = 2H^ mode . From Eq. (|30() in the limit {—krj) <C 1, Ci is related with the scalar field 
perturbation as 

«2J0 = ^:^, (69) 

in both branches. 

In the following, we express the curvature perturbation on the brane as the perturbation of the scalar field in this 
limit in both branches. 

— branch (e = 1) 

In this limit, the = 2H^ mode of 'R-{b){m^=2H^) is, 

'^m^^=^H^)=-^^2nH + ^y (70) 

where we have used the fact that Ji/2i~kri) (— fc77)^/^-y/2/7r as {—kr]) 0. Then we can relate the perturbation of 
the scalar field and the curvature perturbation as 

n,., - ^-Un4i^^, <-e« --^(^|f^, (71) 



which is the same as the standard four-dimensional result (see Eq. JB^TJ), except for the overall normalization of 
the effective gravitational constant. If we take the limit rcH — s- oo, K4_„off coincides with K4 and we recover four- 
dimensional results. We can also show that 'R-^b) = ~^(b)- 

+ branch (e = 1) 

In this hmit, for Hrc ^ 1, the = 2H^ mode of 'R-(b)(m^=2H'^) becomes 

- -^(2i7,,_i) (Hrc-iy 

where we have used the fact that 

a{-kri) 21n(-fc77), (73) 
P{-kri) 2kr] as - kr] ~* 0. (74) 

Strictly speaking, the asymptotic value of a{—kr]) and j3{—kri) contain an arbitrary integration constants (see 
Eqs. HA7|I and IjASp . It was shown in Ref. 0, [sil that this degree of freedom corresponds to a radion and dark 
radiation. In the following, we neglect this contribution. 

Then, for Hrc 7^ 1, we can relate the perturbation of the scalar field and the curvature perturbation as 

On the other hand we get 

A(b) = '^X+c«^H- (76) 
Again, if we take the limit Hrc ^ oo, we recover the four-dimensional results. 
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B. Self-accelerating background 

Similarly, for Hr^ — 1, 'R-{b)(m^=2H^) becomes 

C 

^(b)(m2=2H2) = (-2 + ln(-fc77)), (77) 

where we have used the fact that 

a{-k7j) ^ (ln(-M)2, (78) 
f3{-kri) -2fc?7- 21n(-fc?7)(-fc77) as - fcry ^ 0, (79) 

as well as Eqs. (|75)l and JTJ. In terms of the perturbation of the scalar field, this is given by 

T^ib) = + H-kv))H- (80) 

We also get 

A(b) = 0. (81) 

The solution for TZ^b) is diverging on large scales limit {—krj) — )■ 0. In the next section, we study whether this means 
the physical instability on the brane or not. 

C. Effective theory on a brane 

In order to see the origin of the difference between + branch and — branch, it is useful to consider effective covariant 
gravitational field equations on the brane derived by projecting the five-dimensional Einstein equations and applying 
the Israel junction conditions with reflection symmetry at the brane |36l l37| . 



G,,, = (167rGr,)2n^, - E^,, (82) 



where 



Tf_,u = — (87rG') ^Gf_,u, (83) 
n^. = -\T^.»f,^ + ^r^^f^. + ^ [if^pf'^P - (f„")2] (84) 

and E^i, is the trace-free projection of the five-dimensional Weyl tensor. In the background, we take E^,y = by 
assuming a Minkowski bulk. However, perturbations of E^^ are not necessarily vanishing. For perturbations, the 
effective equations are 

1 \ f . 2fc2 1 
^ 2h7J - ^^W) - ^^0)/ = ^<5pi=, (85) 

2 - ^) (7^(b) - H^(b)) = 4(^5q - l^^CiSqE. (86) 



^^^^ 



Substituting the solution for 'R-(b) and ^(b), we find the perturbations of Weyl tensor vanish in — branch. In + branch, 
5pE = Sqe = 0, but there is a non-trivial anisotropic stress; 

5t^e = -\ ,J"'^ Hr.^l, (88) 

[Hrc - 1) 



^""E = \^ilogi-kT^)-2)S^, Hr, = l. (89) 
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Since a oc {—kr])~^, even though ln(— fc?7) diverges as —krj 0, the anisotropic stress from perturbations of E'^^ does 
not diverge in the hmit —krj — s- 0. 

Since gauge invariant perturbations of Weyl tensor do not diverge, the divergence of the curvature perturbation in 
longitudinal gauge Eq. H8U|I docs not imply the instability of the spacetime itself and the divergence is caused simply 
by a bad choice of the gauge. 

In order to show this, let us consider a four-dimensional gauge transformation, 

r]^r]~e\ x' ^ x' - e\ (90) 
By choosing dri^ = e'', the metric perturbations are transformed as 

ds^ = a2{-(l + 2i(b))d?72 + [(1 + 27^(b))5,J + ^(6) ^,j]dx'dx^, (91) 

where 

= A(b)-a^e''-(9^1na)e^ (92) 
7^(b) = 7^(b)-(9^1na)e^ (93) 
Eib) = -e. (94) 
For Htc 7^ 1, we can eliminate the singular part in TZ^b) in limit Hrc 1 by choosing 

e - '^^^+-ff^Hr,~l)2kH 
Then the resultant metric perturbations are given by 



{-kri)5(j3. (95) 



7^(^,) = = -nX^^^—5(j3, (96) 



d,E,b, = «Uc«(-^'y)(^^^^0. (97) 



For Hvc = 1, we can also eliminate the growing part of Ti-Qj) by choosing 

^^-<50(-fc??) (ln(-A:r;) - 1) . 



Then, the resultant metric perturbations become 



7t(b) = = -K^77<50, (99) 



k|0 

At late times —krj 0, E^b) behaves as 



drAb) - -^(-M(ln(-M-1)^0- (100) 



E^b) « {-krif ln(-fc77) 0. (101) 
Thus we can find a gauge where all perturbations remain small and so there is no gravitational instability. 



VII. CONCLUSION 



In this paper, we studied inflaton perturbations confined to a de Sitter brane with induced gravity in a five- 
dimensional Minkowski spacetime. 

For a vacuum brane, the spin-0 mode appears as a discrete bulk mode with m? = 2H^ in the -I- branch, while in 
the — branch there are no normalizable solutions for the spin-0 modes. Since there is another discrete bulk mode 
(helicity-0 mode of spin-2 perturbation) with mass = H^{3Hrc — l){Hrc)~^ in the -I- branch, there is a resonance 
between the spin-0 mode and the helicity-0 mode of spin-2 perturbation for iJr^ = 1. In this paper we introduced 
inflaton perturbations on a brane. Then an infinite ladder of discrete modes with in? = —2(2^ — 1){£ + 1)H^ and 
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= —2i{2£ + 3)H^ are excited. Since there is a mode with m? = 2H^ regardless of the value of Hrc, in the + 
branch the resonance inevitably appears. We obtained the solutions for the curvature perturbation on the brane Ti-ib) 
and studied their behavior. 

At high energies, Hvc — *■ oo, we have confirmed that, in both branches, the four-dimensional general relativity 
solutions are recovered from the infinite sum of the modes. This results remind us that in the case of massive gravity, 
if we introduce a cosmological constant A(4), the vDVZ discontinuity disappears if the limits to ^ is taken for 
fixed A(4) |3lL l32L Is^. However, this argument holds only for Hrc oo. For a finite Hvc, if we do not fix — fcry and 
take large enough — fcry, we can no longer perform the summation of the infinite ladder of the modes and the theory 
deviate from four-dimensional general relativity. In fact, on small scales, for a fixed Hr^ we can see that the junction 
condition becomes the same as the gravitational field equations in Brans-Dicke theory where the brane bending mode 
acts as the BD scalar. 

Then, in order to identify the BD parameter for the effective BD theory, we compared the solutions for Brans-Dicke 
scalar field and the brane bending mode. We have confirmed that, on sufficiently small scales, gravity is well described 
by the Brans-Dicke theory with uj = —3Hrc in + branch and u; = iHrc in — branch. This is consistent with the 
previous work [TsLllflll and the existence of the ghost in -I- branch is confirmed. 

We also studied the large scales perturbations. In — branch, the solutions agree with four-dimensional GR with 
a modified gravitational constant. In + branch, there is an additional contribution to 1^(1) due to the resonance, 
which diverges on large scales for Htc = 1. We identified this contribution as the effect of the Weyl anisotropic stress 
and showed that the anisotropic stress itself does not diverge. In fact, we can find a suitable gauge where all metric 
perturbations remain small, so the resonance does not lead to a gravitational instability on the brane. 

We make comments on future applications of our results. In — branch, if Vf. <C small scale perturbations 

a/k <C can be described by the four-dimensional Brans-Dicke theory. If the perturbations approaches to Tc, 
the gravity becomes five-dimensional and we expect significant effects from the coupling to five-dimensional metric 
perturbations Unlike the Randall-Sundrum model where the small scales perturbations are always coupled to 
five-dimensional perturbations and a quantum vacuum state is hard to be specified 0, HE we can specify a 
vacuum state without ambiguity based on the four-dimensional BD theory on sufficiently small scales. Then we can 
estimate the effect of the coupling to five-dimensional gravity without ambiguity. 

Finally we should emphasize that our analysis is limited to linear perturbations. It has been shown that we need to 
take into account the non-linearity of the brane bending modes before the gravity becomes non-linear. Actually, we 
can recover four-dimensional general relativity due to this non-linearity 26, 40, 41, 42, 43, 44, 45]. Ref. Ig] suggested 
that this non- linear scale becomes very small at high energies by the factor 1 — e2Hrc, which is the same factor seen 
in Eq. (|54|l . The non- linear interactions of the brane bending mode in de Sitter spacetime deserve further study. 
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APPENDIX A: DERIVATION OF f^„2=2H2 IN + BRANCH 

In order to obtain the mode corresponding to — 2H^, we must prepare the solution other than N{y), since the 
junction conditions cannot be satisfied by this solution. We find that depending on whether Hvc = 1 or not, there 
arc two types of the solutions. 

Hrc^l 

First, we consider the case Hrc 7^ 1. In this case, we assume the following form of the solution: 

^(m^-=2m) = (1 + Hy)P{t) + (1 + Hy) ln(l + Hy)Q{t). (Al) 
By substituting Eq. (jAl|) into the junction condition (|33() . we obtain the following relation; 

(1 - 2Hrc)HQ{t) - rc{2H^ - m^)P{t) = Ci^{-k7^yi Ji{-h^). (A2) 

It is worth noting that in Eq. \A2\ . since P{t) is no longer {—kri)^'^/^Ji{—kr]), m? is not a constant but a differential 
operator. On the other hand, we find that the equation of the motion in the bulk is satisfied if the following conditions 
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hold: 

{2H^ ~m'')P(t) ^ -H^Q{t), (A3) 
{2H^ - m^)Q{t) = 0. (A4) 

By eliminating P{t) from Eqs. (|A2|I and ljA3(l . we get 

H{Hrc - 1) 



and by substituting Eq. (|A5|) into Eq. ljA3|) we obtain 



where 



Pit) = ^'^ 2H{Hr,-l) (-fcry)^ ^"^"^^^) sin(-fc,7) + /3(-M cos(-fc,7)}, (A6) 



- A;?7 

);(-fc?7) = / d{-kff){-kf]y^ s\\\{~2kff), (A7) 



/3(-/c7y) = / d{-kf]){-kf]y^{cos{-2kff) -I}. (A8) 
ifr, = 1 

Next, we consider the case Hrc = 1. In this case, we assume the following form of the solution: 

f^(™^=2//^) = (1 + Hv)P{t) + (1 + Hy) ln(l + Hy)Q{t) + (1 + Hy) (ln(l + Hy)f S{t). (A9) 
By substituting Eq. (jA9|) into the junction condition (|33() . we can obtain the following relation: 

- HQ{t) - l(2i/2 - m')P{t) = Ciy|(-fc77)"^ Ji(-fc^). (AlO) 

It is worth noting that in Eq. IjAlOp . since P{t) is no longer (— fcr;)^"^/^ Ji {—kr]), m? is not a constant but a differential 
operator. On the other hand, we find that the equation of the motion in the bulk is satisfied if the following conditions 
hold: 

{2H^ - m^)P{i) ^ {m^ -3H^)Q{t), (All) 
{2H'^ -m'^)Q{t) ^ -2H^S{t), (A12) 
{2H^ ^m^)S{t) = 0. (A13) 

By eliminating P{t) and Q{t) from Eqs. ifnTjl . (|Il2)l and l|IT3)l we obtain 

S{t) = -^(-M"' sin(-fc,7). (AM) 
By substituting (|A14() into (|A12() . we obtain 

Qit) = ^ (Z^{"("^^) M-kv) + Pi-kv) cos(-fc77)}, (A15) 

where a{—krj) and P{—kri) are defined as Eqs. I|A7|I and ljA8|l . Furthermore, by substituting ljA14|l and (|A15|I into 
(|A11|I . we get 



Ci 1 
2iJ (-fcr;) 

where 

1 



Pit) = -77^r^:^{"(-M sin(-fc7;) + P{~kTj) cos{~kf])}, (A16) 



krj 

a{-kf]) = a{~kr]) + ^ I d{~kf]){~kfiy'^[a{~kfi) sm{~2kf]) + l3{~kf]){cos{~2kf]) + 1}], (A17) 



1 



" — A;/; 

fi{-kf]) ^ p{-kri) + ^ I d{-kfi){-kfi)-^[a{-kfi){cos{-2kfi) -1} ~ (3{~kf])sm{-2kfi)]. (A18) 
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APPENDIX B: METRIC PERTURBATIONS IN 4D BRANS-DICKE GRAVITY 

In this ap pen dix, we briefly summarize the resuhs for cosmological perturbations in four-dimensional Brans-Dicke 
theory (see jig for a review of BD theory) . 

In this theory, we start with the foUowing action: 

-CbD = \/~^ i^R - LU-^g^'^dfj^ipd^ip + imattor^ , (Bl) 



where is a Brans-Dicke scalar field and is a Brans-Dicke coupling constant. From Eq. (|B1|I we can derive the field 
equations, 

2^G^, = T^,+T;f,-2(5^,n- V^V,)(^, (B2) 
Dip = C'T, (B3) 
V^T'^" = 0, (B4) 



where C is a constant related with lu as 



For perturbations, the metric is taken as 



+ (B5) 



ds^ = -(1 + 2A)dt^ + aitfil + 2'R)8ijdx'dx^ . (B6) 
When we consider a canonical scalar field with a potential V{4)) as the matter, the field equations reduce to 

Qu^dip + Aipo[-m{nA - n') - v^n] = -(j)'^A + (j)'S(t)' + v,^a^S(j) - enSip' + 2v^dip, (b?) 

Aipo {T-LA -n') = (j}'5(j} + 25p' - 2nSp, (B8) 

2S(p{2n' + n^) - 4:ipo[{2n' + n^)A + ha' - n" - 2Hn' + ^^^{a + n)] 

= (P'^A - (P'6(f>' + V,^a^S(p - 2,595" - 2HSp' + 2W^Sip, (B9) 
ipo{A + n) = -6^, (BIO) 
dip" + 2ndp' - V^dip = 2C^[(j)'A - (jj'dct) + 2V,^a'^d(l)], (Bll) 

where the prime here denotes the derivative with respect to conformal time ry and Ti — a' /a 

If we require Eqs. (|B7|) . (jB8p . IIB9II and (|B10|) reproduces the results of GR in the limit oi 6p ^ and uj — * 00, the 
value of pq is set to be (po = {1/2k1). For the scale much smaller than the Hubble radius and the scalar field evolves 
slowly in the background, Eqs. (|B7|) . IjBSj) . ljB9|) and IjBlOp can be simplified as 



-AipoV^n = (j)'6(j)' + 2V^Sip, (B12) 

-AipoTZ' ^ (j)'6(j) + 2dip', (B13) 

AipoU" -2poV^{A + n) = ~(j)'5(j)' + 2V'^5lp, (B14) 

ipo{A + n) = -5ip. (B15) 

Then the equation for curvature perturbation is obtained as 

n" - V^n = —V^Sip. (B16) 
2p}Q 

If we take general relativity limit, we get 

■R" - V^n = 0, (B17) 

and its Fourier modes of TZ are expressed as 

7^ = ClCOs(-/^^7) + C2sin(-fcr;), (B18) 
where Ci and C2 are integration constants [l^l . The lapse function A is obtained as 

A = -n= -Ci cos{-kri) - C2 sm{-hj). (B19) 
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If we take a large scale limit, that is, —krj — > 0, from Eq. (|B19|I . we can obtain 

n^C\^ const. (B20) 
From Einstein equations, the curvature perturbation is related with the scalar field fluctuation as 

n = (B21) 
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